Primordial perturbations in a non singular bouncing universe model. 
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We construct a simple non singular cosmological model in which the currently observed expansion 
phase was preceded by a contraction. This is achieved, in the framework of pure general relativity, 
by means of a radiation fluid and a free scalar field having negative energy. We calculate the power 
spectrum of the scalar perturbations that are produced in such a bouncing model and find that, under 
the assumption of initial vacuum state for the quantum field associated with the hydrodynamical 
perturbation, this leads to a spectral index = —1. The matching conditions applying to this 
bouncing model are derived and shown to be different from those in the case of a sharp transition. 
We find that if our bounce transition can be smoothly connected to a slowly contracting phase, then 
the resulting power spectrum will be scale invariant. 



PACS numbers: 98.80.Cq, 98.80.Hw 

I. INTRODUCTION 

For more than two decades, inflation has been the 
only available paradigm to solve the standard cosmologi- 
cal problems of flatness, homogeneity, and monopole ex- 
cess. It also predicts, as a bonus, that primordial fluc- 
tuations, assumed to be of quantum origin, could be en- 
hanced to the level required to trigger large scale struc- 
ture formation, with an almost scale-invariant spectrum. 
To date, no model has ever come close to challenging this 
impressive list of successes. 

Inflation cosmology suffers, however, from a few prob- 
lems of its own, whose seriousness is largely a matter 
of opinion. For instance, in a typical realization, the 
underlying parameters (mass and coupling constants of 
the inflaton fleld) must be assigned "un-natural" values 
in order to reproduce the observed temperature fluctu- 
ations in the Cosmic Microwave Background Radiation 
(CMBR) which the mechanism seeds. However, such a 
fine-tuning can be accounted for in various realistic mod- 
els. 

The inflation paradigm is also endowed with two spe- 
cific problems, conceptually much more serious, that may 
ultimately be related, namely the meaning of the trans- 
Planckian [0 perturbations and the existence of a past 
singularity ra. Concerning the latter, many ideas were 
discussed, among which the Tolman Phoenix universe , 
and many others in the seventies |^ , and recently revived 
under the name "ekpyrotic" ||] in the somewhat different 
context of superstring Q inspired brane cosmology 
This model, however, was the subject of many criticisms , 
both from the string ||] and cosmological points of 
view. In its latest version Wm. moreover, the model also 
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contains a supposedly actual singularity p^ . 

Quantum cosmology, in the framework of the Wheeler- 
de-Witt (WdW) equation, exhibits bouncing solu- 
tions fl^ which can be interpreted as truly avoiding 
the singularity, even for flat (/C = 0) spatial sections, 
a possibility strictly forbidden in classical General Rela- 
tivity (GR), unless |l4| some exotic material, with neg- 
ative energy density violating the Null Energy Condi- 
tion (NEC) is introduced, which most cosmologists 
are reluctant on doing. Such a bouncing universe model 
provides a solution to the horizon problem by geodesi- 
cally completing the manifold in the past, and avoids the 
monopole formation if the bounce takes place at a tem- 
perature below that of Grand Unification (GUT); this 
class of models does not however address the question 
of flatness and one must assume /C = from the outset. 
Moreover, in such a context, the trans-Planckian issue 
simply does not exist because the initial conditions for 
the perturbations can be imposed during a phase where 
the universe is as close to the Minkowski spacetime as 
one wishes, without ever passing through a Planck phase. 
Thus, it could be a natural competitor to the inflationary 
paradigm, and it is therefore of interest to estimate the 
primordial perturbation spectrum that it can produce. 

In a previous work ]l6| , we examined the stability of 
a bouncing universe dominated, at the bounce, by an 
exotic hydrodynamical perfect fluid. We showed, by 
computing scalar perturbations using the gauge invari- 
ant Bardeen potential, that such perturbations grow 
unboundedly either at the bouncing point, or at the time 
when the NEC was violated or restored, thereby contra- 
dicting the hypothesis of low amplitude first order per- 
turbation theory jl^, Such models are thus incom- 
patible with observational data, e.g., CMBR data 
according to which these first order effects indeed still 
dominate (over nonlinear effects) at large scales. 

The next to simple possibility consists on using a scalar 
field instead of an exotic fluid. The purpose of this paper 
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is to exhibit a toy model in which a radiation fluid is cou- 
pled to a negative energy free scalar field that is supposed 
to dominate the universe for a limited amount of time, 
during which the bounce occurs ||2l|] . Our universe thus 
comes from a low-density, radiation dominated, contract- 
ing state, passes through a bounce, and connects again 
back to the usual Hot Big-Bang phase |2^. In this con- 
text, we shall be concerned with the scalar perturbations 
induced during the transition between the collapsing and 
expanding phases. 

In section II we set the various constraints our model 
needs to satisfy, and we explain how it can be made 
phenomcnologically reasonable. Then we calculate, in 
section [II, the power spectrum of the perturbations by 



matching the relevant solutions in the various regions of 
interest, and we compare the results with numerical cal- 
culations. Contrary to what one would naively expect 
from a fluid analysis fl^ , we find that scalar perturba- 
tions are perfectly well behaved all along. 

Setting vacuum initial conditions for the quantized hy- 
drodynamical perturbations deep in the low-density ra- 
diation dominated phase, we find that the relevant spec- 
trum of perturbations, at last horizon crossing during the 
expanding radiation dominated era, has a spectral index 
rig = — 1. It is thus incompatible with observational data. 
As it is a model dependent result, further investigations 
of more realistic models [^3[ , from the point of view 
of particle physics, need to be done p5[ |. 

While other models yield a scale invariant spectrum by 
making use of various assumptions , the present cal- 
culations are made with a specific model where the tran- 
sition through the bounce is made with an exact solution. 
This allows one to obtain, qualitatively and numerically, 
the transitions in the Bardeen potential and its deriva- 
tive through the bounce, yielding indications on what 
kind of matching conditions p7| should be proposed for 
perturbations passing through a general bounce. We ob- 
tain the perhaps not so surprising result that the Bardeen 
potential changes sign through the bounce, even though 
its derivative is continuous, contrary to the case of a 
sharp transition This result will be discussed 

in more detail in section in which we present a way 
to obtain a scale invariant spectrum for the scalar per- 
turbation by connecting our bounce and radiation dom- 
inated model with a slowly contracting phase. If such a 
four-dimensional and singularity-free model can be con- 
structed, it will be able to reproduce all the available 
observational data while avoiding most of the questions 
raised by the inflation solution. 

This paper only deals with the scalar part of the per- 
turbations, which we show does not yield a spectrum 
compatible with the data. The tensor part was al- 
ready calculated in Rcf. where the spectral index 
= rig — I = 2 was obtained, and is therefore unable to 
reproduce the data. Finally, we will not be concerned 
with vectorial (rotational) perturbations even though, 
contrary to the usual inflationary case, one could think 
that those have no reason to be a priori negligible with 



a time symmetric scale factor. However, the Universe is 
torque-free |Q since at least the nucleosynthesis epoch 
that occurred at a redshift of z , ~ 3 x 10^. Hence, 

iiucl ' 

the present relative contribution for the vectorial per- 
turbation, which scales as a^^ [|l8j, is expected to be 
Sy <C 10^^^, independent of the scale k at which it is 
evaluated, and hence observationally irrelevant. 



II. THE MODEL 

We shall consider a very simple toy model for which we 
demand the following conditions to hold. First of all, we 
want general relativity to be valid for all times. We also 
impose that at late times, the model should reproduce 
the standard hot big bang case, i.e. there should exist 
a time in which radiation dominates. This implies in 
particular that we assume some amount of radiation to 
be present in our model. We also restrict our attention 
to the spatially flat situation. Finally, the model should 
have a bouncing phase. This means, given that there is 
already some radiation present, that, in the context of 
GR, there must exist some other fluid having negative 
energy. In particular, for the special case at hand for 
which the spatial curvature /C = 0, this means that the 
Null Energy Condition (NEC) must be violated at some 
time near the bounce p^ . 

Realizing such a model is in principle feasible with just 
another fluid, e.g., some stiff matter with negative en- 
ergy, namely one for which the equation of state reads 
p = p < 0. However, it was recently shown that 
such an approach will lead to an overproduction of large 
inhomogcneities at various different times, breaking the 
cosmological principle hypothesis long before nucleosyn- 
thesis. Such an approach is therefore not applicable, and 
we must resort to the next-to-simple possibility, namely a 
free massless scalar fleld which is known to reproduce the 
stiff matter fluid behavior at the background level 
The action we shall start with thus reads 



S = 



f 



16nG 



-gd^x, (1) 



where R is the curvature scalar, e the energy density of 
the radiation fluid, and 4> the scalar fleld. We assume that 
the background metric takes the standard Friedmann- 
Robertson- Walker form 
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with 7] the conformal time. The cosmic time, t, is then 
obtained as the solution of the equation ad?] — dt once 
the scale factor a{rj) is known. Note that throughout this 
paper, we assume the background curvature to vanish, 
JC — 0. In this context, it has to be a particular choice: 
this category of models does not indeed solve the flatness 
problem. 

Varying the action (|^) with respect to the fluid and 
flelds yields the background dynamical equations 



me = 0, ip" + 2nip' = 0, 



(3) 



where Ti. = a' /a, e and Lp are the background space- 
independent values of the radiation energy density e and 
scalar field (j) respectively, and a prime denotes a differ- 
entiation with respect to the conformal time -q. These 
background equations imply 



i2' 



(4) 



where c and d are constant. The energy density of the 
scalar field is given by 



f_ 
2a2 



c 

2^' 



(5) 



and as such it is dominant when a is small and negligi- 
ble when o is very large. These solutions, together with 
Friedmann equation 
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lead to the bouncing solution 



a(?/) = ao\/ f + 



(6) 



(7) 



where the minimum scale factor oq and the characteristic 
bouncing conformal time 770 solely depend on the relative 
quantities of energy density in radiation and scalar field 
at some given time: Cq = c^/{2d) and rj^ = / {2d^ i"^^) . 
In what follows, these two parameters will be considered 
as the relevant ones. 

Before turning to first order perturbations of this back- 
ground, which is the subject of the following section, we 
want to emphasize a point of stability of this model re- 
lated to the "wrong" sign chosen in Eq. (|l|). Indeed, an 
expansion of Eq. (nl) with respect to 
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with [g , tp] the classical part and (h, S(p) interpreted 
respectively as gravitons and scalar particles in a semi- 
classical approach, will inevitably lead to two different 
kinds of instabilities, each arising at a different order in 
perturbation. The first one, with which we shall deal 
later since it is actually the one responsible for the large 
scale structure formation in this model, is second order 
in perturbation (first order in the equations of motion) 
and goes essentially as oc h'^'^dfj,S(j)d„(f. This term is ab- 
sent in ordinary Minkowski space, but is present in the 
cosmological setup we are considering because of Eq. (^) 
in which the classical part of the scalar field varies with 
time and thus behaves as a source for the production 
of gravitons and scalar particles. As it originates in a 
derivative couphng [see FIG. |l|— (a)], the characteristic 
time scale of this instability is that of the classical scalar 
part, in our case the typical cosmological timescale. 

The second instability that must be discussed is much 
more serious, even though, at first sight, it looks innocu- 
ous because of a higher order in perturbation: it is the 




FIG. 1: Diagrams leading to instabilities in the theory (|l|). 

(a) : the dynamical instability whereby the energy contained 
in the scalar field can be used to produce semi-classical per- 
turbations, later to be identified with primordial fluctuations. 

(b) : Vacuum instability. As this process in non zero, the vac- 
uum can spontaneously decay into a pair of negative energy 
scalar particles and a positive energy graviton. 



same term as before, but with the classical part replaced 
by a first order perturbation, namely oc h^^'^ d^5<j)d^5(l) 
[see FIG. |f]— (&)]• The presence of such a process means 
that the vacuum can spontaneously decay into a pair of 
negative energy scalar particles and a graviton, and, due 
to this fact, the energy levels are not bounded from be- 
low. This sounds like a catastrophe, and even more so 
because the only available timescale comes from the cou- 
pling constant, i.e., the Planck time. However, it is clear 
from the figure that the process probability amplitude A 
is A(x p^/AP, with - ~ IQis GeV the Planck 
mass and p the momentum at the vertex. Such an am- 
plitude therefore becomes important when the charac- 
teristic scale p^^ is comparable to £pj. At this point, it 
should be argued that the model of Eq. (|l|) is understood 
as an effective low energy theory which must be imple- 
mented with a cutoff scale much larger than the Planck 
one: as one reaches the Planck energy scale, the theory 
is expected to break down into a completely different one 
such as, e.g., quantum gravity or supcrstring theory. As 
a result, for cosmological purposes, one can safely ignore 
this instability and concentrate on the production of cos- 
mological perturbations. 



III. LINEAR PERTURBATION SPECTRUM 

In what follows, we shall consider perturbations stem- 
ming from the model (|l]) , making use of the gauge invari- 
ant formalism [|l^ ^ . As there are no anisotropic stress 
perturbations in this model, the most general form of 
metric perturbations on the background given by Eq. (|^) 
reads, in the longitudinal gauge, 



ds^ = a2(ry) [(I + 2^)dTj^ - (I - 2$),5y dxMx^] , (9) 

where $ is the gauge invariant Bardeen potential [ pT| . 
Setting also 
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(f> — ip{r]) + S(l){x,f]) and e ^ e{i]) + 6e{x,f]), (10) 

one obtains the radiation fluid current conservation and 
Klein-Gordon equation respectively in the form 



Se' + mSe = |e (3$' + a-^\/^&) , 



(11) 



where a is the gauge invariant fluid velocity potential, 
and use has been made of the relation 5e = ^Sp between 
the energy density Se and pressure Sp fluctuation. Ein- 
stein equations yield, after a bit of algebra Q, 
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x=r]/77o 
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+ 7^$ = {-ip'5(j> + faa) , 



2n' + n^)^ = ^ii^ 

_ _ (12) 

Simple manipulations of Eqs. (O) and (112) permit 
to eliminate the radiation fluctuation in favor of the 
Bardeen potential through the relation 

<^>i' + 6n^i+ [2 {n' + 2re) + -f^a'Sek, (13) 

where, from now on, we assume a Fourier decomposi- 
tion of each variable A into its components deflned 
through 



d^a 



(27r)3/2 



•"^(|xU), 



(14) 



where A^ only depends on the amplitude k of the wave- 
vector k. 

The dynamical equations for the Bardeen potential and 
the fluctuations of the scalar field therefore decouple from 
the radiation fluid perturbations and are then expressible 
solely in terms of themselves as (making use of the back- 
ground Einstein equations) 

<^l! + mH + \k'<i>k^-el^ip'54>i, (15) 



and 



b{! + 2U54>L + fc'(50fc= 4^'$^. (16) 



We shall now investigate the solution of these equations 
in order to get the perturbation spectrum such a bounc- 
ing model predicts. 



A. The relevant phases in the perturbations 
evolution 

In order to investigate Eqs. (^5|) and ([l6|), let us write 
them in terms of the variables Wfc = a?'^k E^nd Wk = 0(50^ 



FIG. 2: Potentials for the parametric oscillator equations giv- 
ing the dynamics of the Bardeen potential and the scalar field 
perturbations. The full line shows the potential for the vari- 
able associated with [see Eq. (^)], the dashed line the 
potential for the scalar field Wk [Eq. (|2^)], and the dotted 
line, showing the value of (riok)^, indicates visually the dif- 
ferent regions where the different approximations hold. The 
points xi and X2 are the matching points for these two fields. 



Using Eq. (Q), they read 
1 



;P-2{n' + 2T-e) 



Uk 



and 



{e -H^- H') Wk = 4— {u'^ - Huk) 



Hwk) , 
(17) 

• (18) 



Each of these equations can be seen as an inhomogeneous 
equation, i.e., one with a source term not depending on 
the function itself. Asymptotically, since the scale factor 
grows like jr/j, it can be checked explicitly, by means of 
an expansion in powers of 77"^ for Uk and w^, that the 
source terms in both Eqs. ( |l7| ) and ( |l^) are small for^ 
r] ^ rjQ. Therefore, in this limit, Eqs. (|15D and (|l6|) can 
be reduced to the usual parametric oscillator equations 
for the variables Uk and Wk, namely 
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and 



a 
a 



Uk = 0, 



Wk = 0. 



(19) 



(20) 



The potentials {a^)"/a^ = 2{n' + 271^) = 2/(772 + 772) for 
Uk and a" /a = H' + = rj^ / {if + 'Ha)'^ for Wk are shown 



^ This is valid not only when ^ 00 but also when < 6{W + 
2?i^) as long as 770. We shall return to this point in more 
detail below. 
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on FIG. g on which we also define the variable x = r]/r]o, 
as well as the various corresponding matching values Xi 
and X2- 

We shall be interested in the cosmologically relevant 
limit fc ^ 1. However, as \ri\ — > oo, the fc-dependent 
terms in Eqs. ( |l9| ) and ( ]20|) become important. More 
precisely, when k\r]\ > v6, i.e., when \x\ > xi = 
VQ/{kr]o) ^ 1, the solutions of the above equations in 
terms of $fcand 5(/)/jCan be written in terms of Hankcl 
functions, namely pi|. 



-3/2 



1] 



-1/2 



*(l)^3/2M)^ 



^(2)^3/2^ 



(21) 

where lu = fc/\/3. 

When the potential terms dominate over the fc- 
dependent terms, which for Uk is the case so long as 
fc|r/| ^ -\/6, or 1 ^ \x\ <^ xi, and for Vk when k\r]\ <^ 
\/fc??o, or 1 ^ <C X2, where X2 = lU/^qo, the ze- 
roth order solutions for Eqs. (19) and ( PQ ) read, when 
xi <C a; <C a;2, 
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FIG. 3: The bounce functions fi as functions of x = rj/rjo- 



with A, B, C and D arbitrary constants. The bounce 
functions fi{x) are found to be 



/i(^) 



(l + a;2)2^ 



f2{x) 



l-x^ 



2(l + a;2)2^ 



(26) 



/4 
a-2dr, + 0(fc77) « Ai-A2^" 



(22) 



I' 



M?/ + 0(fcr;) « Bi 



i?2^. (23) 



In fact, solution (|2^) needs amelioration because around 
I a; I = X2 the source term of Eq . ([iq ) become important. 
This is not the case of solution (|22|) because for 1 ^ \x\ < 
X2, the source terms of Eq. (|l5|) are still negligible, even 
taking into account the corrections of (p3|). Fortunately, 
for what follows in the next subsection, only the solution 
( ^ ) will be needed. 

Near the bounce, when the potentials and c/a are of 
order l/r/^ and ao/(?7o^pi)i the source terms become im- 
portant2, but the terms proportional to fc2 are still neg- 
ligible. In this situation, one can neglect altogether the 
fc2 term in Eqs. (^5|) and (|l6|), yielding the solutions 



and 



<^ '"'"'^ =A + Bh{x)+Cf2{x), (24) 



,50 =D + Bh{x) + CUix), (25) 



^ The constraints on the values of ag and rjg in order for the bounce 
to happen at a scale much larger than the Planck scale and a long 
time before nucleosynthesis takes place are £pj ^ ao'?0 ^ 10* cm. 
Hence, if one chooses rjo of order one, the potentials are of order 
one, and c/a 3> 1 near the bounce. 



V2 

'{1 + x^y 



/4(X) 



y2(l + a;2)2 



(27) 



and are displayed on FIG. ||. 

These solutions will be used to match the asymptotic 
solutions through the bounce. 



B. Matching the solutions, and the power 
spectrum 

In the limit i] ±oo {<=i^ a — > ±oo), i.e., very far 
from the bounce, the Universe is radiation dominated, so 
that the coupling term in the left hand side of Eqs. ( [l7| ) 
and (|l^ can be neglected, as it was explained in the last 
subsection. From Eq.(|2l]), and for 77 — — oo, the Bardeen 
potential and the scalar field perturbation respectively 
scale as l/a2 and 1/a. From Eq. ( p^ one can see that 
the fluid perturbation Sek goes like 1/a*. Therefore, from 
Eq. (12), we find that the scalar field and its perturba- 
tion are irrelevant, in this regime, for the evolution of the 
Bardeen potential with respect to the radiation fiuid. For 
this reason, one can conclude that the appropriate quan- 
tum gauge invariant variable to be used must be the same 
as the one defined in Ref. for the quantum treatment 
of hydrodynamical fluids perturbation theory, which, in 
the case of pure radiation, is related to $ by. 



2^2 



(28) 



where P = H'^ - H' and z = a^/3P/H. Similarly, the 
gauge invariant quantum variable connected to the scalar 
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field perturbation given in Ref. |19| | is given by 

uJk = a[5<j)+{ip/ny^]^Wk. (29) 

It is interesting to note that the quantum field v leaves 
the oscillatory regime at the same conformal time as S(j) 
does, and that neither of them do so at horizon crossing. 
This is a peculiarity of our model due to the fact that at 
the time at which the quantum fields leave the oscillatory 
regime, the space is almost radiation dominated, but not 
quite. 

Imposing the initial vacuum state for these quantum 
variables implies that we can set 



3l/4g-ifc(j,-j,i)/\/3 



Vk = 



and 



2k 



-iMv-Vi) 



Wk 



2k 



at rj —^ — cxD, with rji and rjj two a priori arbitrary confor- 
mal times, having no influence on the subsequent evolu- 
tion. From the solutions (|^) and these initial conditions, 
one can write the Bardeen potential (f>feand the scalar field 
perturbation J^^at ki] <^ — \/6 (or |a;| 3> a:i 3> 1) as 



2aof]^k^/^ 



(^:^ + i}jc-^Hv-m)/V3^ (30) 



and 



aorjV^k 



(31) 



We are interested in calculating the power spectrum 



- 1,3 



k'\<^k 



(32) 



evaluated at the time when <i>fc returns to its oscillatory 
regime, i.e., at cc = cci. As we shall see later, the values of 
6(f>k 'm the different phases of perturbation evolution are 
not necessary to calculate $feat x = xi. Hence, we will 
forget about (5(/)fefrom now on. 

Looking at Eq. (pT|), one can see that the first matching 
must be imposed when k'^/3 = 2(H' -I- 27^^) = {a'^Y'/a'^, 
for Uk = cL^^k As k is very small, this happens when 
I??! ^ 1 (where we can ignore the source terms). Match- 
ing the solution (pw with solution (22) at the point 



kr] « —^/6 (or x ~ —xi) yields 



y/k 



35/42^2 



i{V2+krii)/V3 



(33) 



and 



(1 - 372z)e'(^+^''^)/^. (34) 



2?73fc5/2 



A. 



The solution (|2^) is valid up to the point where x is of 
order one, when we approach the bounce. Differentiating 
Eq. (15) twice and making use of Eq. ( p^ ) as well as 
the background equations, we obtain the following fourth 
order equation 



(IV) ,,, 



+ 20 {H' + 2H^ 



+ GHk^^L + -k^ [fc2 + 4 (W + 2%^)] $fe= 0, 



r 



(35) 



whose analysis indicates that the solution (|22| ) is valid 
for \x <C a;2, up to |a;| « 1, where the source term in 
Eqs. (1^) become relevant but k is still completely neg- 
ligible. We have to match this solution with the other 
relevant solution in this region, namely, the bounce solu- 
tion ([2^), which have three arbitrary constants. As con- 
firmed by the numerical analysis below, we have to set 
C = because the function it multiplies goes like l/x^, 
which should appear, and dominate, in Eq. (^) in a re- 
gion where both solutions are valid. Hence, we only need 
to determine B and A. That is why we do not need to 
calculate the evolution of (50^ in order to determine such 
constants; this also explains why the initial conditions we 
assume for the scalar field perturbation are irrelevant for 
the final power spectrum. 



The solution (gj) will propagate the Bardeen potential 
to the other side of the bounce, to the region where x is of 
order one. As we are in a region where k is negligible, the 
point of matching will be chosen tobea; = — A^^— 1, 
where N does not depend on k but is large'^. 



^ We consider A'^ large but not large eno 
order in the expansion of /i in Eq. 



I to neglect terms of 
) . If we neglect such 
terms, we loose the effect of the bounce in the evolution of the 
perturbations. Also, considering N = 1, without approxima- 
tions, would not change our qualitative results and the power 
spectrum. 
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The result of the matching reads 



and 



B 



8m 



A2 



3a, 



(36) 



(37) 







At a; = iV ^ 1, on the other side of the bounce, 
these solutions must be matched with a solution simi- 
lar to Eq. (E3), namely 



— Ci — C2 



3a4?73 ' 



yielding 



Ci=Ai 



(38) 



(39) 



term in Eq. ( 38 ) is the constant Ci : as we are now back 
to a regular expanding universe, the other term is a de- 
caying mode which rapidly becomes negligible. In C'l, 
the dominant term when /c <C 1 is the one proportional 
to A2 which goes as k~^/'^ , while the other is proportional 
to ^/k. Hence, we get 



krjf) (and correspondingly u) = Lorjo), and write the system 
as 



da;^ 



Ax d^k 
x-^ + 1 dx 
2x dXk 



1 da; 



k^Xk = 



V2 dXk 
:2 + 1 dx 
AV2 d$fc 
+ 1 dx 



(42) 

relations in which Xk = ^pjiJ^fc (recall that x = rj/rjo), 
subject to initial conditions, far in the limit x — > —00, 
given by Eqs. ( ^ and (^) with r]i = r]j = , namely 



33/4a /V3 



2x'^k^/'^ V kx 



ikx / y/S 



(43) 



and C2 = A2. For the power spectrum, the important and 



Xk 



-ikx 



'2k 



(44) 



fc3|$fe(-xi)pcxP|v42| 



(40) 



yielding a spectral index rig = — 1. 

One can then define a transfer function between "Hori- 
zon exit" and "Horizon re-entry" as the ratio of the power 
spectra at the corresponding two different times. It is 
given in the case of our bounce by the relation 



fc^|$..(-xi)p 6 



(41) 



This transfer function essentially depends on the behav- 
ior of the scale factor at both times, as well as on the 
nature of the bounce itself. It is represented on Fig. ^ 

Let us now check all these approximations through a 
numerical examination of Eqs. ( |l5| ) and (p^). 



Numerical calculations 



The system ( [I5[]l6| ) can be solved numerically for any 
value of k. For that purpose, we also include the char- 
acteristic conformal timescale 779 in the wavenumber k — 



where we have defined the only free dimensionless param- 
eter a = ^piy^/ag. In all the figures, this parameter has 
been arbitrarily fixed to the value a = 10~^; the conclu- 
sions do not, however, depend on this value, which acts 
as a simple normalization constant. 

The solution of Eqs. (^ ) for the square of the Bardeen 
potential |$fcp is shown on the bottom panel of FIG. ^ 
for various values of the wavenumber, renormalized with 
the bounce characteristic conformal timescale, fc, rang- 
ing from 10~® to ~ 1 on the figure as a function of the 
renormalized conformal time y = kx = krj. All calcula- 
tions are started far in the radiation dominated epoch, 
for y — —100, where the boundary conditions hold. This 
is verified as, indeed, for small enough values of fcjyo, 
behaves as 77^"', as expected. It can be checked 
that, as discussed above, in the long wavelength limit, 
the Bardeen potential starts with a negligible constant 
part and a growing, oc 77"^, mode, for — xi ^ x <C 1, 
which then connects to the /i part while crossing the 
bounce, and then connects back to the usual growing 
and decaying modes, although the new constant part has 
now acquired a piece from both modes of the previous 
epoch. 



Once the system ( |4^ ) is solved, one can easily compute the value of the Bardeen potential at horizon crossing, 
namely for x ~ 1/^, i.e., 77 ^ 1/fc, or y ~ 1. This provides the spectrum shown on the top panel of FIG. ^. It is clear 
on that figure that for small values of fc, the behavior of the power spectrum is indeed a power law, which we checked 
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indeed corresponds to rig = —1. Also shown is a comparison between various cases of interest, namely the vacuum 
case for which the initial conditions given by Eqs. ( ^3|) and ( ^ ) hold, the gravitational vacuum case for which Eq. ( |4^ ) 
still holds, but with Eq. (^4|) replaced by X^r' = 0, and finally the decoupled case for which the coupling between 
<I>fe and Xk is made to vanish, i.e. for which the left-hand side of Eqs. (|42| ) is arbitrarily set to zero. The curves 
corresponding to either vacuum or gravitational vacuum initial conditions are seen to be almost undistinguishable, 
showing that, as expected and discussed in the previous section, the final spectrum for the gravitational potential does 
not depend on the initial conditions for the scalar field perturbations. The decoupled curve shows that, for fc <C 1, 
if one were to neglect the bounce duration and apply some matching conditions by brute force, one gets the same 
spectral index — —1, but with a normalization that is wrong by many orders of magnitude. The situation is even 
worse for intermediate scales for which even the index is wrong. 



Vacuum initial conditions 
Decoupled case 
Sip -0 




10" 



kr\=10' 



-40 







20 



40 



FIG. 4: Top panel: spectrum of scalar perturbations, i.e., fc as function of the wavenumber fc, normalized with 770 as 

indicated. The long wavelength part of the spectrum, as expected, is well fitted by a power-law with spectral index — — 1. 
The full line is for vacuum initial condition for Sip, the dotted line is with Sip = Sip' = at the initial time (j/ini = krjini = —100 
in the numerical calculation), and the dashed curve represents the fully decoupled situation for which ip' is assumed negligible 
all along. Bottom panel; Time evolution of the gravitational potential |$fe|^ for different wavelengths. 



On FIG. ^ is shown an enhancement of the region sur- 
rounding the bounce itself. This figure shows that the 
real and imaginary parts of both the Bardeen potential 
and the scalar field perturbation connect, respectively, 
with the bounce functions /i and /s, thereby confirming 
the prediction (7 = 0. 

The fact that the Bardeen potential only connects to 
the odd bounce function fi{x) suggests that in the limit 
in which the bounce duration 770 can be neglected, one 
may apply the following junction conditions across the 
surface at which the bounce is located (77 = 0) 



= =0, 



(45) 



where [A]^ = A{—r]g) — A{+n^) is the jump in the geo- 
metric quantity A. On FIG. are shown the time evo- 
lution across the bounce of the quantities involved in 
Eq. (^) . For a fixed surface thickness ij^ , and in a way 
independent of this thickness, the relevant quantities are 



indeed conserved and can be safely used. This is of course 
true only in the particular example presented here, but it 
can also be conjectured to apply for a symmetric bounce 
in general. In fact, taking the ekpyrotic model Q with 
our matching conditions (45), one obtains a scale invari- 
ant spectrum. 



IV. CONCLUSIONS 

We have presented a cosmological model in which a 
bounce takes place in the framework of pure general rel- 
ativity. This is achieved by assuming that, at some stage 
after a contracting phase, a negative energy free scalar 
field became important. Performing a bounce with such a 
scalar field, instead of an ordinary hydrodynamical fluid, 
permits to regularize the perturbation, which otherwise 
grow unbounded near the bounce p6|. We derived the 



9 



T(k) 

-- k'm 



FIG. 5: Transfer function for the bounce model. Full line: ra- 
tio of the squared gravitational potential amplitude between 
horizon exit and re-entry. The dashed line shows the same 
multiplied by to emphasize the power law behavior ob- 
tained in Eq. (^). 
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FIG. 6; Real (left) and imaginary (right) parts of the Bardeen 
potential (top) and the scalar field perturbation (bottom) as 
a function of a; = v/Vo for k = 10~*. 



last horizon crossing spectrum and obtained, both ana- 
lytically and numerically, a spectral index — —1 in the 
long wavelength limit, therefore ruling out such a model 
as a competitor to the inflationary paradigm. However, 
our study of a concrete bouncing model allowed us to 
obtain some intuition on what happens with perturba- 
tions when they pass through a bounce. First of all, 
the bounce acts indeed as a "pump field" for perturba- 
tions. Secondly, the field which produces the bounce in 
the background solution, and its perturbations, is not rel- 
evant for the evolution of the Bardeen potential in almost 
the whole history of the model, except near the bounce 
itself, where it becomes very important for the power 
spectrum amplitude, although not for the spectral in- 
dex. Finally, usual matching conditions are not valid for 



-2x10 
-2xl0' 
-2xl0' 
-3xl0' 




2 3 4 5 
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2xl0'^ 
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FIG. 7: Real (left) and imaginary (right) parts of the geomet- 
rical junction variables 7i$fe(top) and Tl$k+ ^L, (bottom) as 
a function oi x — ri/r)o for k — 10~*. 



transitions through a bounce. In fact, even for the back- 
ground metric, such conditions are not valid since the 
Hubble parameter H = Ti/a changes sign through the 
bounce, by definition. Through our bounce, the Bardeen 
potential also changes sign, and what happens to be con- 
tinuous is the combination 7i$. Inspired by our concrete 
model, we suggested matching conditions to be applied to 
general models where the bounce is no t specified, which 

Of course, these 



are spelled out in the subsection |III C , 
suggestions must be checked within other concrete exam- 
ples, or through a more general formal analysis. 

The model we have discussed is admittedly over sim- 
plistic. We may be confident in its latest part describ- 
ing the radiation dominated epoch, which we know have 
taken place in our Universe, and accept that it may pro- 
vide a reasonable description of an immediately preced- 
ing bouncing phase. There is, however, no reason to be- 
lieve, even in the case of a bounce, that the evolution of 
the Universe should have been symmetric in time. One 
can instead set up a contracting phase with a different 
scale factor, assuming at some stage some form of entropy 
production, to end up with enough radiation before the 
bounce, yielding a scale factor of the form (^ that ul- 
timately connects back to standard cosmology. Such a 
model should originate in a realistic underlying particle 
physics theory. 

Let us briefly discuss an example, which is reminis- 
cent of the ekpyrotic proposal 1^, with a few differences. 
First, the model we have in mind would be purely four 
dimensional and does not intend to address the flatness 
problem. Second, such a model would be effectively sin- 
gularity free. Finally, we would not need to impose arbi- 
trary p^ ] matching conditions across the bounce to ob- 
tain the required mixing in the growing and decaying 
modes before and after the bounce, since its specific form 
would be known. 

More precisely, a model satisfying the abovementioned 



10 




a slowly contracting phase, a oc {—ijy, with < p 1, 
connected to the phase examined in the previous sections, 
as the one shown on FIG. ||. When the perturbation 
in the Bardeen potential crosses the horizon for the first 
time, the dependence rj^^ in Eq. (^0[), stemming from the 
fact that the universe is supposed to be radiation domi- 
nated at that time, would be substituted by a dependence 
77~^P, i.e., almost independent of k when kjj ^ 1. Doing 
calculations along the lines of those presented in Sec. Ill, 
the scale invariant spectrum follows. 
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FIG. 8: Scale factor for connecting a slowly contracting phase 
to our model. If such a four dimensional model was effectively 
constructed, it would produce a scale invariant spectrum of 
perturbations and would thus become a promising competitor 
to more usual inflationary models. 
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